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The problem to be addressed in this work 
is that of modeling nondimensional force 
and moment aerodynamic coefficents over 
the entire subsonic envelope. The particular 
application discussed here is the Z force 
coefficient for the F-18 High Angle of Attack 
Research Vehicle (HARV). 
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System Identification for Aircraft 


A priori knowledge: 
Rigid Body Dynamics 
Aerodynamics 



Wind Tunnel Data 
Flight Test Data 


High Fidelity 
Mathematical 
Description 
of the 

Aircraft Dynamics 


In general for high fidelity math models 
of aircraft dynamics, some a priori knowledge 
is combined with experimental data. The 
experimental data generally supplies the 
information on the aerodynamics particular 
to the aircraft under study. 
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Aircraft Equations of Motion 


Fore# 

Equations 


q-SC, + T 
U = - QW + RV + gx + — m 


V = -RU + PW + gy + 


TS Cy 
m 


w = - pv + qu +g I + 3 4r i - 


Shown here are the three nonlinear 
equations of motion for translational motion 
of an aircraft. Moment equations are 
omitted for simplicity, but the same 
discussion can be applied to them as well. 

The objective is to find a model for the 
aerodynamic coefficients, like Cz, 
in terms of the aircraft states and 
outputs and their derivatives, and the controls. 
This would allow prediction of dynamic 
behavior when such models are substituted 
into the equations above, along with a model 
for the thrust. 
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Aerodynamic Modeling 

C.-C.Ca, 

where 5 = Sp- . “ = ## 



C t = Co + Cj a + C2 M + C3 (J + C4d+C5^ + C«8n + C75f 
+ Cg CtM + C9 Ct* Q + ... ? 

Kev Questions; 

• What terms should be Included ? 

• What are the parameter values ? 



♦ When Is a given expansion valid ? 


Part of the a priori knowledge is the 
functional dependence shown on the 
first line of this slide. Using this information, 
a polynomial model is postulated which is 
analogous to a truncated multivariable 
Maclaurin series expansion for Cz. This 
approach gives rise to the three questions 
posed at the bottom of the slide. 
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Aerodynamic Coefficient Expan sion 


C z = c 0 + cia + c 2 M +c 3 d+C4a+c 5 ^ + c«5 11 + C76f 

+ eg a M + C9 <x 2 Q 


i iQiiallv difficult to accur ately estimate parameters 
In this expansion f rom experimental b$C9USg l 


• Experimental data does not contain sufficient Information 

• Control system or the experiment correlates the regressors 


• The model 


Assuming now that the number of terms to 
be included has been determined somehow 
(perhaps using stepwise regression), and the 
range of validity of this expansion has also 
been determined, there remain further 
difficulties in accurately estimating the 
parameter values (the c's) in the model. These 
difficulties are listed on this slide. 
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( j Aerodynamic Coefficient Expansion ^ 
in Terms of Orthogonal Functions 




where in general 

• & = ( a, M, Q, 5, $i, 5„, 5f ) for i = 1, 2, ..., n 


• §| T = 0 for i * j 



the I th term in the model 
^ depends only on c z and 5, 



If instead of ordinary polynomials in the 
expansion for Cz, orthogonal polynomials with 
the properties shown on this slide were used, 
the difficulties in accurately estimating the 
model parameters disappear. This is because 
the computation of the value of the model 
parameter for term i depends only on the 
measured data (Cz) and orthogonal function i. 
This follows from the expansion 
for Cz in terms of orthogonal polynomials 
and the properties of orthogonal polynomials. 
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Model Structure Determination 


• Minimum predicted square error Is the criterion for which 
orthogonal functions should be included in the model : 




The orthogonal functions to 
be included in the model can be 
determined in a straightforward way 
using the predicted square error (PSE) 
criterion, where n is the number of 
terms in the model and N is the number 
of data points. Since each term in the 
orthogonal function expansion is 
independent of all the others, the decision 
of whether or not to include each individual 
orthogonal function term can be made 
based on whether or not its inclusion 
reduces PSE. This can be done for each term 
sequentially and without regard to any 
other terms already in the model. 


395 





Finding Ordinary Polynomial Models 

It is possi ble to ; 

• Generate orthogonal functions based on the data using 
a technique similar to Gram- Schmidt orthogonalizatlon 


• Determine the orthogonal function expansion by 
decrease in fit error vs. additional terms 


* Decompose each orthogonal function without ambiguity 
In terms of ordinary polynomials, e.g. 





* Arrive at an ordinary polynomial model with 

adequate structure and accurately estimated parameters 



Details of how to generate multivariate 
orthogonal functions based on the data 
can be found in the reference given on the 
first slide. In that reference, it is also shown 
that each orthogonal function can be 
decomposed precisely into an expansion in 
terms of ordinary polynomials. Once the model 
structure determination is done in terms of 
orthogonal functions, the expansions for 
each included orthogonal function in terms 
of ordinary polynomials is substituted and 
common terms are combined, resulting finally 
in an ordinary polynomial model with 
adequate model structure and accurately 
estimated model parameters 
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Full Envelop e Expansion 
for the F-18 HARV 
from Wind Tunnel Data 



c z 


ci + c 2 a + c 3 o 2 + C4 6 s + c 5 a 3 + c 6 afc + c 7 a 2 k , + c 8 M 
+ ( c 9 + c i0 a + c n M + c 12 a 2 + c«aM + cu&M 
+ cis a 3 + ci6« 4 + ci 7 cPM + ci*a 4 M ) Q 

+ (ci9 + c 20 a + c 2 iM+c 22 rf 

+ c 22 aM + c 2 4a 2 M + C25M 2 + c 2 6d 3 ) « 

+ (C 27 + C 2 gO ) 5n + (C29 + C30<* + C31M ) 8f 



Valid fori 


- 4.0 deg £ a £ 90 deg 
0.2 £ M £ 0.9 

- 24.0 deg £ & £ 10.5 deg 

- 3.0 deg £ 8» £ 33.0 deg 

- M dee < 8f S 45 i 0deg 



This is the result of using the 
technique developed here to model the z 
aerodynamic force coefficient for the 
F-18 HARV over the entire subsonic 
flight envelope, based on a wind tunnel 
database. 
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A prediction case for the polynomial model 
is shown here, using independent variable 
time histories taken from flight data. The 
prediction is excellent, with slight error in the 
time period from 13 to 23 seconds. This was 
traced to the fact that the wind tunnel database 
zeroed certain components of the 
Z force coefficient when angle of attack 
exceeded 40 degrees, because of lack of data. 
The polynomial model, on the other 
hand, extrapolated reasonably, as shown on 
the lower right plot. 
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Features of the Model 




• Largo tables of aerodynamic data reduced to 
a small number of parameters 

• Smooth aerodynamic functions and partial derivatives 

• Better Insight on the aerodynamic functional dependence 

• Potential for automated simulator updates 

• One model for the entire flight envelope 

J 


Along with the improved accuracy 
of model parameters and excellent 
prediction capability inherent in 
the modeling approach discussed here, 
this slide outlines the some features 
of nonlinear aerodynamic models 
generated with this technique. 
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Conclusion 


> \ 


Accomplishment^ 

• Developed and demonstrated a full envelope 
nonlinear aerodynamic modeling technique 

[Applications | 

• Flight Simulators 

• Global Stability and Control Analysis 

• Dynamic Analysis 

• High Angle of Attack Control Design 

J 


The fact that the aerodynamic model 
for a full subsonic envelope can be 
made compact and smoothly 
differentiable while retaining excellent 
predictive capability has potential 
utility in a wide range of activities. 
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